study of scaling in the fractional quantum Hall effect 
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In the composite fermion model of the fractional quantum Hall effect, composite fermions ex- 
perience, in addition to the usual potential disorder, also a magnetic flux disorder. Motivated by 
this, we investigate the localization properties of a single fermion in two dimensions, moving in the 
presence of both potential and magnetic flux disorders, but with a non-zero average magnetic field. 
It is found that the exponent characterizing the divergence of the localization length is not changed 
upon the addition of the flux disorder, provided it is not too large, suggesting that the transitions 
between fractionally quantized Hall plateaus belong to the same universality class as those between 
the integrally quantized ones. 

PACS numbers: 73.40.Hm 



In the absence of a magnetic field, an arbitrary weak 
random potential localizes all states of an electron in two 
dimensions [1]. In contrast, existence of extended states, 
separated by bands of localized states, is crucial for the 
explanation of the quantum Hall effect [2]. Numerical 
studies [3,4] have made a convincing case that the local- 
ization length of an electron, confined to two dimensions 
and exposed to a strong magnetic field, diverges at a sin- 
gle energy near the center of each Landau level, and the 
divergence is characterized by an exponent a: 

^=\E-E,\-", (1) 

where a 7/3, independent of the Landau level (LL) 
index. No analytical derivation of the value of the expo- 
nent yet exists. Experimentally, the width of the tran- 
sition region (e.g., the width of the conductance peak) 
in the integer quantum Hall regime (IQHE) regime has 
been found to vanish in a range of temperature as , 
where /3 0.42, independent of the LL in question [5]. 
This is consistent with the theoretical value of a in a vari- 
able range hopping theory [6], which obtains /3 = 1/a (in 
other words, writing /3 = l/za, this work shows that the 
dynamical exponent z is equal to unity for the Coulom- 
b interaction). Direct measurements of the localization 
length are also in agreement with a 7/3 [7]. 

At first glance, the nature of localization would ap- 
pear to be much more complex in the fractional quan- 
tum Hall effect (FQHE) regime, since interactions are 
an intimate part of the physics here and cannot be ne- 
glected. However, in the composite fermion model [8], 
the essential part of the interelectron interaction is in- 
corporated through the formation of composite fermion- 
s, which themselves can be treated as non-interacting to 
a good approximation. An effective single-particle de- 
scription of the strongly correlated electron state thus 
becomes possible in terms of composite fermions. This 
paper studies the scaling of the localization length in the 
FQHE regime in the framework of a non-interacting com- 



posite fermion model. While the dynamics of composite 
fermions is similar to that of electrons in a weaker mag- 
netic field, an essential difference is that the composite 
fermions also experience a random fiux disorder in addi- 
tion to the usual random potential disorder [9-11]. The 
main conclusion of the paper is that the exponent a is 
the same for non-interacting composite fermions as for 
non-interacting electrons, provided the amplitude of the 
random fiux disorder is smaller than the average mag- 
netic field. This is consistent, following Ref. [6], with the 
experimental result that the exponent for the transition 
from 1/3 to 2/5 is the same as for the transitions between 
IQHE states [12]. 

The result is not entirely unexpected, and earlier works 
have also anticipated that the FQHE transitions belong 
to the same universality class as the IQHE transitions. 
Jain, Kivelson and Trivedi [13] argued that the compos- 
ite fermion (CF) wave functions describe the physics even 
in the presence of disorder; this leads to a law of corre- 
sponding states relating the FQHE states to the IQHE 
states, and makes it plausible that the localization length 
of composite fermions is the same as that of electrons to 
which their wave function is related, leading to the same 
exponent. The validity of the CF wave functions in the 
presence of disorder was demonstrated only in some lim- 
ited cases, however. The present work investigates the 
problem from a different point of view, and sheds further 
light on the issue. 

According to the CF ansaiz, in a range of filling fac- 
tors in the lowest LL, electrons find it energetically fa- 
vorable to capture an even number of vortices of the 
many particle wave function. The bound state of an 
electron and vortices behaves as a particle, called the 
composite fermion. The vortices produce phases as the 
composite fermions move around, which partly cancel 
the Aharonov-Bohm phases originating from the external 
magnetic field, and, as a result, the composite fermions 
experience an effective magnetic field given by 
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B* = B - 2mp<l)o 



(2) 



where B is the external field, 2m is the number of vortices 
carried by composite fermions, (fiQ = hc/e is the magnetic 
fiux quantum, and p is the electron (or CF) density. In 
this manner, the problem of strongly correlated electrons 
at B is mapped on to the problem of weakly interact- 
ing composite fermions at B* and the FQHE transitions 
of electrons are mapped onto the IQHE transitions of 
composite fermions. One might be tempted immediately 
to conclude that the same exponent a should charac- 
terize the divergence of the localization length. Howev- 
er, there is one important difference between the IQHE 
transitions of electrons and composite fermions: while 
for electrons the magnetic field is constant, the effective 
magnetic field for the composite fermions fiuctuates in 
space, since a potential disorder produces density varia- 
tions, which, through the above equation, leads to fiuctu- 
ations in the magnetic field experienced by the composite 
fermion [9,10]. Motivated by this picture, we consider a 
single composite fermion in an average magnetic field B* , 
in the presence of a random potential as well as a (static) 
random fiux disorder. 

In the following, we briefiy outline our model and tech- 
nique. Our numerical sample consists of a very long two- 
dimensional strip of a finite width (M) square lattice with 
nearest neighbor hopping. Periodic boundary condition 
in the width direction is used to avoid complications aris- 
ing from the edge states. The disorder potential is mod- 
eled by an on-site white-noise potential Vim (*' denotes the 
column index, m denotes the chain index) ranging from 
— W/2 to W/2. The magnetic field appears through the 
complex phase of the hopping term. The strength of the 
magnetic field is characterized by the fiux per plaquette 
((/)) measured in unit of (fiQ- It is the sum of an average 
field (fia and a random field distributed uniformly between 
—4>rl2 to -\-4>rl2. The Hamiltonian of this system can be 
written as: 
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For a specific energy E, a transfer matrix Ti can be easily 
set up mapping the wavefunction amplitudes at columns 
i — I and i to those in columns i and i + I, i.e. 
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Hi-E -I 
I 



(5) 



where Hi is the Hamiltonian for the ith column and I is 
an M X M unit matrix. Using a standard iteration algo- 
rithm [4], we can calculate the Lyapunov exponents for 
the transfer matrix Ti. The finite-size localization length 
^m(E) at energy E for finite width M is then given by 
the inverse of the smallest Lyapunov exponent. In our 
numerical calculation, we choose the sample length to 
be over 10'* so that the self-averaging effect automati- 
cally takes care of the ensemble statistical fiuctuations. 
We then obtain the thermodynamic localization length ^ 
using the standard one-parameter finite-size scaling anal- 
ysis [14], according to which the renormalized finite-size 
localization length Xm I M can be expressed in terms of 
a universal function of M/^, i.e.. 
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where f(x) cx l/x in the thermodynamic limit (M co) 
for localized states while approaching a constant (~ 1) 
when ^ diverges. The critical exponent a, given in Table 
I, is obtained from the best fit of ^ = \E — Ec\~" . The 
value of Ec is obtained using the method of Ref. [15]. 
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where < im;jn > indicates nearest neighbors on the lat- 
tice. The amplitude of the hopping term is chosen as the 
unit of energy. The gauge is conveniently chosen so that 
the inter-column hopping does not carry a complex phase 
factor {i.e., tim,i+i,m = —1) and the effect of magnetic 
field shows up only through the phase factor of the intra- 
column (inter-chain) hopping term. If the random fiux 
in a plaquette cornered by (im), (i + 1, m), (i -|- 1, m -|- 1) 
and (i, m -|- 1) is (fiim, then 
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FIG. 1. Thermodynamic localization length { for average 
field (l>a = 1/7 and zero random field (l>r = with disorder 
W = 1.0 (circles), W = 2.0 (squares) and PF = 3.0 (triangles). 



We focus below on the localization behavior in the low- 
est Landau level, i.e. on the 1^2 transition of compos- 
ite fermions, which corresponds to the 1/3 2/5 transi- 
tion of electrons. Throughout this work, the average field 
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is fixed at = 1/7. Fig. 1 presents tlie results for the 
disorder dependence of the thermodynamic localization 
length ^ with zero random field. The results are clear- 
ly consistent with the notion that all states are localized 
except at a single energy, Eg, which is determined to be 
Ec ~ —3.2 from the method of Ref. [15]. We note that Eg 
can also be obtained from numerical fitting, which gives 
almost identical answer. The main effect of the disorder 
is to change the localization length; interestingly, in the 
present problem, increasing W causes the ^ to increase 
in contrast to the usual intuition. [16] 



10000.0 




0.1 I ' ' ' ' ' 1 

-4.0 -3.5 -3.0 -2.5 

E 




0.1 I ' • ' • ' 1 

-4.0 -3.5 -3.0 -2.5 

E 

FIG. 2. Thermodynamic localization length { for average 
field (l>a = 1/7. (a) W = and random fields (l>r = 0.1 
(circles), (l>r = 0.2 (squares) and (l>r = 0.3 (triangles); (b) 
W = 2 and random fields (l>r = (circles), (l>r = 0.1 (squares), 
(l>r = 0.2 (triangles) and (l>r = 0.3 (diamonds). 

Next we compute the localization length ^ as a func- 
tion of the random magnetic field (for several values of 
disorder), shown in Fig. 2. For a given energy, ^ increas- 
es for increasing random field (fir, indicating that the role 
of the random field is qualitatively similar to that of the 
random disorder. The critical exponent a is again calcu- 



lated by fitting £,{E) <x \E - Ec\~", and listed in Table 
I below. Two exponents are given for each value of W 
and (fir, one obtained from the lower energy side and the 
other from the higher energy side. The error shown in 
Table I for each critical exponent represents the error of 
the least square fit for ^ versus energy plot. Another con- 
tribution to the uncertainty, of order ~ 0.3, comes from 
the various possible choices of points close to Eg used for 
the fit. The exponents are consistent with 7/3 within our 
numerical uncertainty, demonstrating that the addition 
of random magnetic field does not alter the universality 
of the transition. 
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FIG. 3. Scaling function for average field (l>a = 1/7 and 
disorder W = 1.0 with random field (l>r = (circles) and 
(l>r = 0.6 (squares). 

Until now, we have restricted our study to the situa- 
tion where the amplitude of the random field is smaller 
than the average field ((fr = 0.3 being a borderline case, 
with (/)r/2 only slightly larger than (fa = 1/7). Under this 
condition, the chirality of the particle does not change in 
the sense that the cyclotron motion is everywhere in the 
same direction. Once the amplitude of the random field 
exceeds the average field, the cyclotron motion can be 
clockwise or counter-clockwise in various regions of the 
two dimensional system. Fig. 3 plots the scaling function, 
i.e., Xm/M versus M/^ for (f)^ = 1/7 and W = 1 with 
(fr = (circles) and (fr = 0.6 (squares). The curve with 
circles is a typical representation of the scaling function 
for all the cases listed in Table L In all those cases, we 
found the single-parameter scaling hypothesis of Eq.(6) 
is satisfied in the sense that all data fall into one smooth- 
ly universal curve similar to the circles in Fig. 3. On the 
other hand, the data (squares) for (fr = 0.6 clearly do 
not fall into one smooth curve, indicating a lack of scal- 
ing behavior. In particular, there are no data points for 
M/^ 0, suggesting that ^ does not diverges. In Fig. 4 
we plot ^ for W = 1 (circles), W = 2 (squares) and 
W = 3 (triangles). The result is qualitatively different 
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from that in Fig. 2, especially on the high energy side of 
Ec, where the CF eigenstates are much more extended, 
and the localization length fluctuates strongly as a func- 
tion of the energy or disorder. Further study is needed 
to clarify several issues here. 
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FIG. 4. Thermodynamic localization length { for average 
field (l>a = 1/7 and random field (l>r = 0.6 with disorder 
W = 1.0 (circles), W = 2.0 (squares) and W = 3.0 (tri- 
angles). 

Before closing, we discuss an important assumption in 
the above model, namely that of the neglect of the residu- 
al interaction between composite fermions. In addition to 
potential and static flux disorders, the inter-CF interac- 
tion will also produce dynamical flux disorder. However, 
various theoretical and experimental studies have shown 
that the interaction between the composite fermions is 
weak, and can be included perturbatively. Studies treat- 
ing the Coulomb interaction perturbatively in the IQHE 
regime have found that it does not affect the critical ex- 
ponent [17]. This suggests that the inter-CF interaction 
is also unlikely to alter the critical behavior. 
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TABLE I. The critical exponent(a) for several values of 
random magnetic field {4>r) and on-site disorder (IT). 
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